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TRANSIENT RESPONSE OF LOW FREQUENCY YERTICAL ANTENNAS 
TO HIGH ALTITUDE NUCLEAR ELECTROMAGNETIC PULSE (EMP) 


P. R. Barnes and D. B. Nelson 


ABSTRACT 


Several antenna models are used to calculate the tren- 
sient response of top-loaded LF vertical antennas to EMP 
from high altitude nuclear detwnations. Examples are given 
of the short circuit current and voltage across an idealized 
load for a 360 meter high, top-loaded antenna illuminated 
by two representative EMP fields. Goodagreement is obtained 
among all the models used. 


The effect of the angle of incidence 6 is quite pro- 
nouiced. The early time current and voltage at 30° is almost 
twice what it is at 90°, confirming the sin-' 6 behavior. 
On the other hand, end effects occur sooner at 30°, reducing 
the low frequency content. Thus the late time behavior is 
largest at 90°. Ingeneral itwill be true for long antennas 
that effects depending on initial rate of rise will be maxi- 
mized for small 09, whereas effects that dependon total surge 
energy are maximized for broadside incidenc~. 


1. INTRODUCTION 


1.1 Purpose and Content 


This report has been prepared for two purposes: first, to exhibit 
methods for calculating the response to transient electromagnetic waves 
of long thin antennas used for low frequency (LF) radio transmission; and 
second, to give examples of the surge currents and voltages induced in 
such antennas and their coupling circuits by the electromagnetic puise 
(EMP) from high altitude mclear detonations. 

Because LF monopoles are normaily several hundred meters tall, they 
are electrically leng for much of the EMP spectruz- Accordingly, many 
antenna resonances must Ye included to compute the response accurately. 
Because the models used describe the antenna in the frequency domain, 
with time domain answers given via an inverse Laplace or Fourier trans- 
form, a wide range of frequencies is required. To assure that correct 
results are obtained we have employed several theoretical models valid 
in different regions and have compared results in the overlapping 
regions of validity. The simplest of these is the semi-infinite cylin- 
drical antenna erected perpendicular to an infinite ground plane; this model 
gives valid results until the first end reflection reaches the point of 
interest. For later tines the finite length antenna must be considered. 
Two different models are used in this regime. The first of these is a 
classical theory of the monopole, using the work of King” for the 
electrically short antenna and of wu for the long antenna. (We have 
also used a modification of this due to Shen, Wu, and King. ’°) The 
second is the singularity expansion method, developed by Baun® and 
applied to the cylindrical dipole and monopole antennas by Barnes.! 
This application of the singularity expansion method is presently 
developed for broadside wave incidence (electric field parallel to the 
antenna) but it is easy to compute and gives good results for late time. 
The classical method is the more powerful, in principle giving results 
for all time and all angles of incidence. In practice, its reliance on 
a brute force Fourier transform along the real frequency axis makes it 
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Geometry of EMP-Antenna Interaction. 
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prone to many numerical difficulties. We have obtained good results 
from this model by careful selection of frequency points using the 
results from the other models. 

For the epplications we have considered a moropole approximtely 
360 meters tall and 0.8 meters in radius. A top loading array has been 
included, albeit somewhat empirically. The watching network, inductive 
and resistive, resonates with the antenna st 34.5 kHz. Two representa- 
tive EMP fields were considered; each of these was taken at two different 
angles of incidence. The response functions chosen were the short cir- 
cuit current and the voltage across the load. 


1-2 Electromagnetic Environment 


References exist which describe EMP from high altitude nuclear 
8,9 


detonations. For our purposes we idealize the EMP to 4 plane wave 
with the propagation vector in an arbitrary direction. For simplicity 
the electric field is taken to lie in the plane formed by the antenna 
and the propagation vector. This polarization results in a horizontal 
magnetic field. Figure 1.1 illustrates the geometry of the wave- 
antenna-ground system. The angle of incidence is 9A; broadside incidence 
corresponds to 90°. For the examples, we have chosen 9 to be 90° and 30°. 

Tne time dependence of the electric field is idealized to a double 
exponential 


E(t) = E,(eo* - e Pt) . (1.1) 


Two different sets of parameters have been chosen: a short pulse with 
Eo,a, and B given by 


E. = 1.61 x 10® v/m 
a = 5-00 x 10” sec} (1.2) 
B =1.76 x 10° sec™? 


“The term "monopole" is used in this paper to describe a half-dipole 
erected perpendicular to a ground plane. 
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and a long pulse with 


i} 
So 
| 


= 4.20 x 10*V/m 
a = 4.00 x 10® sec-? (3 
4.76 x 108 sec™} . 


D 
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These units correspond to time in seconds and electric field in volts/ 
meter. 


2. ANTENNA ANALYSIS 


2-1 The Low-Frequency Radiator 


The low-frequency (Li) vertical radiator generally consists of a 
vertical tower erected above a radial ground system. Since the physical 
height of the tower is normally less than one-eighth of the wavelength 
of the transmitting frequency, the LF vertical radiator is usually an 
electrically short antenna. The characteristics of an electrically 
short vertical radiator may be improved by adding a top-loading umbrella. 
The top loading tends to increase the antenna's effective height and 
decrease the voltage between the base of the tower and the ground system. 
The geometry of a typical LF verticai radiator is shown in Fig. 2.1. 

When the LF transmitter a’:tenna is viewed as a receiving device, its 
load is the transmitter and coupling circuitry. The antenna is normlly 
connected to the coupling circuit by a short wire-above-ground transmis- 
sion line. The inductive reactance of the coupling circuit and the 
traucmitter is adjusted to cancel the antenna's capacitive reactance at 
the transmitter carrier frequency. Also, the resistance of the coupling 
circuit and transmitter is designed to match the radiation resistance of 
the vertical radiator for maximum power transfer. The value of that 
resistance is normally on the order of 1 ohm. A typical LF transmitter 
output circuit is shown in Fig. 2.2a. 2, is the impedance of the trans- 
mitter and the remaining coupling circuitry. 

When an antenna is considered as a receiving device it is useful to 
employ a Thevenin or a Norton equivalent circuit. The Thevenin and Norton 
equivalent circuits are shown in Figs. 2.cb and 2.2c, respectively. The 


equivalent circuit parameters are defined below: 


V voltage across the load 


I - current through the load 


N 
i] 


load impedance 


Y. - load admittance 
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Fig. 2-1. Typical LF Vertical Radiator. 
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(0) SIDE VIEW-VERTICAL RADIATOR WITH TOP-LOADING UMBRELLA. 
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(5) THEVENIN EQUIVALENT CIRCUIT 


(e) NORTON EQUIV..ENT CRCUIT 


Fig. 2.2. Equivalent Circuits for Determining the LF Antcnna Rcsponse. 
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Z, ~ impedance of the antenna 
oe - admittance of the antenna 


- equivalent voltage source of the antenna and the 
electromagnetic environment 


- equivalent current source of the antenna and the 
electromagnetic enviroment. 


Once the equivalent circuit parameters have been determined, the 
load response can be calculated by circuit theory. The following cir- 
sait relationships apply: 


V I 

V= a. = Pty (2.1) 
a z a L 
LY. V 

I = Y ~ Y = Z = (2.2) 
alee ier ad 

T= Vv’. . (2-3) 


Also, from Fig. 2.2a the load impedance ic given by 


hye 


2,(@) = Jel, * 2, + Jeb, + jal, (2-4) 


For the purpose of calculating the equivalent circuit parameters, a 
thin cylindrical antenna can be used. A thin cylindrical antenra with 
@ noncircular cross section such as a solid tower behaves like a circular 
cylindrical antenna with an equivalent radius. For frequencies of 
interest, the tower behaves like a solid structure. The equivalent 
radius of a tower with an equilateral triangular cross section is*° 

a, = -l2a, (2-5) 


Where 4, is the radius of the circle which circumscribes the equi lateral 
triangle. 


=li- 


2-2 Early Time Model 


The vertical tower radiator can be idealized by an early time zodel 
consisting of semi-infinitely long, perfectly conducting, circular 
cylinder erected perpendicular to a perfectly conducting ground rlane 
as shown in Fig. 2.3. The time domain solution of the infinite antenna 
response is the same 93 that of the finite antenns for early tines, 
i.e., times before tae refiections from the to; of the finite structure 
have effect. 

The Korton equivalent current source for a semi-infinite, perfectly 
conducting, cylindrical structure of radius “a” erected above a perfectly 
conducting ground plane is twice that of an infinite cylindricel antenna 
and is given by” 


ha cE(s) 


T,(s) = C sin@ sk, (es sina/c) (2.6) 


where s is the Laplace transform variable, E(s) is the Laplace transform 
of the incident wave time history, ¢€ is the free space wave impedance 
approximately equal to 120 x ohms, c is the free space speed of light and 
6 is the angle of incidence measured from the axis of the antenna. 

The admittance of the semi-infinite cylindrical antenna separated 
from the ground plane by a gap distance "b” can be derived from that of 
an infinite cylindrical antenna as! 


where 


u = Ve? + (as/c)? (2.8) 


and 
‘= b/a . (2.9) 


K, (u) is the modified Bessel function of the second kind and order i. 
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Fig. 2.3. Semi-Infinite Cylindrical Antenna and the Incident Plane Wave. 


A siaple early tine asyuptotic fora of the voltage response across 
an inductive load can be derived for a double exponential incident wave 
using Equations (2-1), (2-6), and (2-7):4 


v(t) ~ = ae (2.10) 


where the tine history of the incident wave is given by 


E(t) = E, (eo -e pt ) (2,11) 
and time t equal sero corresponds to the tine that the incident vave 
reaches the equivalent driving point of the antenna. From Bq, (2.10). it 
is evident that the first term in the asyaptotic expansion of the voltage 
response across an inductive load is not a function of the inductance. 
This is a reasonable result since at early times (high frequencies) the 
inductance has a near infinite impedance and the first term of the solu- 
tion for the load voltage is equal to that of the open circuit voltage 
response. 


2-3 Moderate and Late ‘ime Models 


For moderate and late tine, the length of the structure and the top- 
loading umbrella have important effects on the response of the vertical 
radiator. The finite length produces end reflections in the time domain. 
The top-loading umbrella adds more resonances corregponding to the 
umbrella-vertical antenna-ground system. It affects the quasi-static (low 
frequency) antenna parameters by increasing both the static antenna 
capacitance and the effective length. 

The finite length thin monopole (or dipole) has been studied exten- 
sively. The methods we use, the singularity expansion method and the 
theories of King, Wu and Shen, are outlined in Sections 2.4 and 2.5. 
However, a rigorous treataent of the umbrella is not extant for the wide 
frequency range required. An empirical approach is used here, matching 
the quasi-static behavior of the radiator. In this approximation the low 
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frequency admittance of the antenna is modified to agree with the quasi- 
static capacitance of the monopole plus umbrella, while a fictitious 
height h is chosen equal to twice the low-frequency limit of the effec- 
tive height of the antenna. 

The quasi-static response of the LF vertical radiator is a first 
order approximation for the moderate and late time response of the 
structure. This can be inferred from the time domain solution by the 
singularity expansion method (sem).! For moderate and late times, only 
the singularities near the origin in the complex frequency plane contri- 
Bute appreciably to the solution. Thus, it is the low-frequency 
singularities that are important for the moderate and late time response. 

Because the modi“ications due to the umbrella differ slightly in 
the two theories used for moderate and late time behavior, they are 
discussed separately in the next two sections. 


204 Application of the Singularity Expansion Method 


A new technique to calculate the dominant electric cylindrical 
dipole and monopole antenna responses to a transient electromagnetic 
plane wave pulse has been described in a previous Pepores! Basically, 
it involves applying the singularity expansion method formulated by 
Baum. This method expands the frequency domain solution in terms of 
its singularities. The time domain solution is the inverse Laplace 
transform of each term in the frequency domain singularity expansion. 

The Norton equivalent circuit parameters for a cylindrical monopole 


antenna can be expressed in terms of their singularities as! 


1, (8) = >> An(5*7n) * Pay 


(8°71 (84%)? + 0,3 +o 2 (2.12) 


M 
Y,(s) = De Me (2.13) 
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The constants ‘nm and o. are the natural damping constants and the 
natural resonant frequencies respectively and are functions of the 
antenna parameters only. The constants A, and Bo determine the ampli- 
tude of the current and are functions of both the antenna parameters and 
the angle of incidence of the transient electromagnetic pulse. E(s) is 
the Iaplace transform of the incident field time history and h is the 
physical height of the antenna. The constants Ka determine the amli- 
tude of the admittance and the constant M is the number of singularities 
considered in the solution. The values of the constant quantities 
defined above are listed in Ref. 7 for broadside incidence and a wide 
range of antenna length-to-radius ratios. They have not yet been con- 
sidered by Barnes for other angles. 

The method employed in Ref. 7 to derive the singularity expansion 
quantities of the monopole admittance is accurate for frequencies equal 
to or less than the third resonance of the antenna. With M = 3 an 
accurate solution can be obtained for late times, i.e., t =>6h/c. If 
five singularities are used, the solution is within a few percent of the 
exact solution for times t >h/c. However, for t < h/c large errors 
could result due to the inaccuracy of the admittance function. 

The modifications necessary to account for the top-loading umbrella 
are indicated in Fig. 2.4. The actual height is replaced by a ficti- 
tious height h, equal to twice the low frequency limit of the effective 
height of the antenna. The equivalent admittance network of the quasi- 
static model is shown in Fig. 2-4b. The admittance of the cylindrical 
monopole hes been modified by the addition of C.» the capacitance due to 
the top-loading umbrella, to that part of the network consisting of R,, 
L, and C, which is associated with the first resonance of the antenna. 

Y is the admittance of the remaining part of the equivalent antenna 


a-1 
admittance network. The total antenna admittance can be written as 


s(C,+C,) 
b ee 6 


+ 2-14 
a a-1 1+ s*L, (C\+ C,) + sR, (s)(C,+C ( ) 


> 


and, of course, the antenna impedance is given by 


ee te 
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(7) PHYSICAL MODEL FOR CALCULATING THE EQUIVALENT NORTON 
CURRENT AND THEVENIN VOLTAGE SOURCES — A CYLINDRICAL 
MONOPOLE OF LENGTH A. 

R, (s) 


Ly 


C, | C, 


(5) EQUIVALENT ADMITTANCE NETWORK OF THE TOP-LOADED 
ANTENNA. 


Yo-4 


Fig. 2.4. Quasi-Static Model of the LF Vertical Radiator. 


<n eed TRAE ee NT ee PET ON NNT Tp ep Re I ee cm ee 


s172 


Zz, ° 5 (2.15) 


ror t = 10h/c. the solutio: for broadside incidence can be accurately 
characterized by only one singularity. If only one singularity is used, 


the Norton circuit parameters are given by 


[A(sty) + Bw, | 


I. = -Sisk(s) 2 2.16 
= (sty)? + uF (2-16) 
and 
esk 
+a (2.17) 
(stv)? + w? soe 


where the "1" subscripts have been dropped for A,, B,, K,, and 7,- 
The current through an inductive load L due to the first singularity 


is 


ShsF(s) [A(s+7) + Be? 
lo SC; (2.18) 


(sty)= + wy + 2s*Kwqh 


The value of K can be computed approximately from the ar.tenra 


static capacitance C,, which includes the top-lcadirg capacitar.ce, by 


Co (7? + a2) 
K oa oa ; (2-19) 


Tor the electrically short antenna 
wilC, >>1l - (2.20) 


ar.d we can also observe from Ref. 7 that for thi:. cyli:drical artenras 


Bo>A (2.21) 


Aarnd 


ea a CaaS 1 a 
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Substituting E4s. (2.19) and (2.20), (2.21) and (2.22) inte Eq. (2.18) 
and performing the inverse Laplace transform gives an approximate late time 


current response to an exponentially decaying incident pulse 


{m’\2 _y? 
i(t) = She, rn Be” * cos w’t (2.23) 
where 
1 

oS (2.24) 

V I, 

yu’? 

LT ee (2.25) 

(7? + wi) 


and time history of the incident wave is 


E(t) = E,eo* (2.26) 


where a >>(7’+W‘). 


The late time voltage across the inductor is given approximately 


by 
=¥'t 
: -8hE,L(w’)* Be 
V= LS ~ gg sinw’t  . (2.27) 


2-5 The Classical Antenna Theories 


The classical theories of the monopole stem from the work of King 
3 of the elec- 
trically long antenna, with later simplications.‘?° The Thevenin 
equivaient of the antenna, shown in Fig. 2.1b, is determined by two 


on the electrically short antenna and the treatment by Wu 


functions of radian frequency w (jw = s, the Laplace transform variable 
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used before): the antenna impedance Z,(), and the effective height 
h, (w) ; Z,(#) is used directly in Thevenin equivalent, and h, (w) deter- 
mines the voltage source Vv, (@) through the forma 


Vv, (#) = ~2h, (@)E, (w) ; (2.28) 


Here E(w) is the component of the free space incident electric field 
parallel to the antenna. Ground reflection (equivalently, the antenna | 
image) accounts for the factor of two in the forma. 

Different functional forms for Z, and h, are used for electrically 
short and electrically long antennas. With k defined by k = ©/c, where 
c is the speed of light, the dividing iine is kh = 1. For kh <1 the 
antenna is electrically short, and the relevant formas are!” 


IW, cos kh 
Z(») =- 7 [sin kh + T(kh) (1 - cos kh) a ie (2.29) 


h (w) — SOS (kh cos 9) - cos kh 2.30 
ei” k sin§® sin kh os 


with 
C = free space impedance = 1201x 
W = radian frequency 
hi 
Va ; 2Zme - 2 
h = antenna height 
&@ = antenna radius 


2 2 : 
2 n2-=)kh - j=(kh)* 
2 n (2h/a) - 3 


For kh >1 the antenna is electrically long and the formlas of wwe are 
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used. Since Wu assumed a time dependerce exp(-iat), instead of exp(j%t), 
the complex conjugate (*) of some of his quantities is required for 
agreement with the rest of this report. Ze and h, are 


23k # m1 5 
Z,(a) = [= (s+ ay] (2.31) 


jnZ, (@) sin 9 


h, (@) kG , - &sin A” 


I 


jc" {kh(2-cos Q)]j ‘ sin [kh(l+cos 9) y (2.32) 
1 - cos 6 ar eoe a 


if - cos [kh(1-cos 6)] . 1 - cos [kh(l+cos si|| 
+ —_ + AE REE TT A I EE cE PL EE NE 
2 1 - cos @ 1 + cos @ 


Here 


2T-T’) sin kh - (2S-S’) cos kh 
T’ cos kh + S sin kh 


U 


ad (A, az Ag) 
yz 0.57722 


ds & 1.6449 


or (2/ka) - 7 


& 
i 


Og = Ng- m2 


O, = 0, + in/2 
N= 226 + mM2kh + 7 - jx /2 
Qp'= Qe + ne 


od + IX Pind 1 1 
si i Pearce Aare] 


A, = (3/2 2) + 7" (a3? - 237) 
As = n(ag/aZ) + 37’ (ag? ~ O57?) 


As = ah 23H (9,7? - 057) 


S = 5 (-A, +A, + As) 

81 = 3 (A, +8 + AS) 

T 3 (-A,- A, +Ag) 

Th = 3 (Ay -AZtAS) 


The formulas of Wu are quite complicated. More recently, two papers 
by King, Wu and Shen *?? developed a simple theory of the long transmit- 
ting and receiving linear antenna which uses mich more intuitive techniques 
than the Wiener-Hopf method of Wu- The resulting formas for 2, and h . 
are correspondingly simpler. A further advantage of this theory is that 
its starting point is the current induced by a plane wave in an infinite 
circular cylinder, our Eq. (2.6). In undertaking the present study we 
were particularly interested in the variation of antenna response with 
angle of incidence 9- Equation (2.6) shows this to be approximately as 
sin7+6. (The apparent singularity at 8 = 0 is removed when finite 
antenna conductivity is included. The infinite conductivity model gives 
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approximately the same results as the finite conductivity uodel with rea- 
sonable conductivity values until 9 is equal to or less than some small 
angle.) ‘This fact means that the standard treatment, broadside inci- 
dence, produces the least current in an infinite wire. As we are in 
general interested in worst case response, it is important to consider 
angles of incidence other than broadside when studying long antennas 
such as the LF monopole. The formulas of Ref. 5 allow easier understend- 
ing of the variation in effective height with 9 than do those of Wu. 

This simple theory leads to computed values of Z,, (w) and h, () which 
are virtually identical to those of Wu. The formulas are 


Ze. = 5197? (2.33) 
* Tx, 
be = F sine T_(0) ov 


with 


Ip = 1,(0) + 2¢,1_(h) 


Jkz ang 
10 ah - ere a] 
a2) C = 2c, + On(kz + V (kz)? + e727) + 7 +5 3/2 


C, = ~onka - 7 


y ~ 0.5772 
C ~ 120% 


C4 = “RI (h)/f1 +RI (2h) ] 
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lt, = Ven) * (Cy, + Cyg) 1 (h) 


tae j exp (jkz cos 9) 
se C sin @ {C, ~ log [(sin 9)/2] + jx/2} 


: RR, (x-6) I, (hb) 1 (2h) - 1, .(-h) R,(6) 
= 1 = R?I_?(2h) 


BR, (8) 1,_(-h) 1 (2h) - 1, (hn) R, (x- @) 
Cso = 


1- R*I_2(2h) 
R= (2c + jx) 


R.(6) =R- a log [(1 - cos @)/2] . 


All of the formulas given above are for a simple monopole above an 
infinitely conducting ground. The actual LF antenna uses a top-loading 
umbrella which affects both Zz, and h,- A rigorous theory of this struc- 
ture would be quite difficult to derive. Empirically it can be added 
to the classical theory in a manner very similar to that used in the 
singularity expansion method of Section 2.4. The change in effective 
height is simulated by a change in the real height h so that h = 2h, (0). 
The change in Z, cannot be made exactly as in Section 2.4, because the 
values of L, and C, associated with the first antenna resonance are not 
readily obtainable from the classical theories. Instead the umbrella is 
treated quasi-statically as a LC resonant structure in parallel with the 
antenna reactance but in series with the antenna resistance. That is, 

Z a is separated into real and imaginary parts: 


ys? 7 


bed 
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Zo) = R, (@) + 5X, (@) (2.35) 
and the reactance of the umbrella is denoted hy X,: 


X,(@) = -1fec, +eL, - (2.36) 


Then the modified impedance Ze is 


xX x 
re at 
Z, = R, + IX +X 

a t 


(2-37) 


Cc, is the known static capacitance of the umbrella, and since the 


horizontel length of the umbrelia is approximately the same as the 


antenra height, lL, is chosen so the umbrella resonates at the first 


antenna resonance: 


2 
Dy, = as (2-38) 


While this treatment of the umbrella is certainly not rigorous, it 
is correct in both the low frequency and high frequency limits and pre- 
sumably mirrors the gross behavior of the monopole-umbrella system. 
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3- ANTENNA RESPONSE CALCULATIONS 


3-1 Introduction 


The short-circuit current and load voltage response of an example 
LF transmitting antenna system to the two double exponential EMP's of 
Eqs. (1-1)-(1-3) are considered in this section. A 1200-foot top-loaded 
vertical tower radiator with an effective radius of 0.75 mis modeled by 
a cylindrical monopole antenna. One effect of the top-loading umbrella 
is to increase the effective length of the antenna from about 366 m to 
450 m. Also, the static antenna capacitance is increased from about 4 nf 
to 14 nf by the top-loading umbrella. Thus the additional capacitance 
due to the top-loading, Cy» is 10 rf. 

The models and caiculational techriques described in Section 2 have 
been employed to calculate the response curves presented here. Plots 
calcuiated by the different methods are presented on the same graph for 
comparison. 


3-2 The Short-Circuit Current 


The short-circuit current responses of the vertical radiator to 
the long and short incident plane wave pulses are shown in Figs. 3.1, 
3-2 and 3.3 with 6 as a parameter. In Figs. 3.1 and 3.2, the early 
time and early-to-moderate time responses are shown respectively. The 
solid curves have been calculated by the classical antenna method and 
the dashed curves have been calculated by the exact solution of a semi- 
infinite cylindrical antenna with Y=1.0. (¥ is defined in Eq. 2.9.) 
The relative error of each dashed curve is less than 0.1%. In Fig. 3.2, 
the dashed curves are plotted for only the time intervals during which 
they are valid for the 1200-foot top-loaded antenna. 

Tue classical antenna method can introduce errors in the early 
time results due to the limited frequency range feasible and numeri- 
cal Fourier inversion used in the calculations. The highest frequency 
considered was w = 0.8 x 10°, which corresponds to a resolution of 


t 12 nsec. However, as can be seen from Figs. 3.1 and 3.2, there is 
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Fig- 3.1. 
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Early Time Short-Circuit Current Response. 
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Fig. 3-2- Early to Moderate Time Short-Circuit Current Response. 
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good agreement between the results calculated by the classical antenna 
method and the semi-infinite antenna method ever for times less than 
10 ns. Thus, the early time errors introduced by the classical antenna 
method are reasonably small. 

In Fig. 3.3, the moderate time short-circuit current responses of 
the LF vertical radiator are shown. The dashed curves have been cal- 
culated for 9 = 90° by the singularity expansion method employing five 
singularities. The solid curves have been computed by the classical 
antenna method with 0 as a parameter. As can be observed from Fig. 
3-3, there is good agreement between the two methods for times greater 
than 1.4 wsec. For eariier times, not enough singuiarities nave been 
used to accurately compute the dashed curves. 


3-3 The Load Voltage 


The voltage responses across a typical LF load connected to the 
1200-foot top-loaded vertical radiator are shown in rigs. 3-4 through 
3-7- The load network shown in Fig. 2.2a with L, = 1.21 mh, L, = 0.2 
mh and 2 = 190 was used forthe calculations. In all the figures the 
solid curves are the results as computed by the classical antenna 
method. The dashed curves are the results as calculated by the semi- 
infinite antenna method in Figs. 3.4 and 3-5 and by the singularity 
expansion method in Fig. 3-7- As can be seen from the plots, the 
solid and dashed curves are generally in good agreement. 

In Fig. 3-6, only the classical antenna method was used to compute 
the load voltage responses. The singularity expansion method normally 
gives accurate results for moderate times; however, the formulation pre- 
sently used by Bernes employs an approximation of the antenna admittance 
which gives inaccurate results for early-to-moderate times. 

In Fig. 3-7, the late time load voltage responses are shown. The 
dashed curves were computed using only one singularity and are reason- 
ably accurate for times greater than 20 usec. As can be seen, there is 
good agreement between solid curves for 8 = 90° and the dashed curves. 


Some numerical errors resulted from the Fourier inversion used by the 
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Fig- 3-3. Moderate Time Short-Circuit Current Response. 
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Fig. 3-4. Early Time Load Voltage Response. 
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Fig. 3.6. Moderate Time Load Voltage Response. 
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Fig. 3-7- Late Time Load Voltage Response. 
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classical antenna method for times greater than 70 usec. However, these 
errors were easily removed by averaging the high frequency oscillations 


introduced by the inverse Fourier program. 
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4. CONCLUSIONS 


4.1 Comparison Among the Antenna Models 


Our main reason for using a number of antenna models was to provide 
sufficient checks on the computed results. The very good agreement 
obtained in all cases indicates that there are no significant numerical 
errors rvesent. Some adjustments were required. For example, to obtain 
the early time results of Figs. 3.1 and 3.4 the classical method required 
frequency domain values up to kh = 120. This includes approximately 75 
antenna resonances and over 700 frequency points for the inverse Fourier 
transform. 

The models describe the LF vertical radiator imperfectly, particu- 
larly in their treatment of the umbrella. Although the early and late 
time results should be quite accurate, the moderate time response of 
Figs. 3-3 and 3.6 could be in error since the models do not include the 
separate resonant effects associated with the physical lengths of the 
tower and the umbrella alone. Such error would likely be in the 
time behavior, the peak values of current and voltage are probabiy 
fairly accurate even for moderate time. The oscillatory wave associated 
with the umbrella is poorly coupled to the LF antenna’s terminals and 
thus would have a small effect. The resonant effect of the physical tower 
height could decrease the first zero crossing time of the response by as 
much as 20 percent. However, the peak of the current and voltage will not. 
be affected since the time of the first peak is sufficiently less than the 
first zero crossing time. Other deficiencies in the models include the 
assumptions of perfect conductivity for antenna and ground plane, treat- 
ment of load elements by lumped parameters, and neglect in all but the 
infinite antenna model of the gap between tower base and ground. Com- 
parison of results shows this last error to be insignificant for the 
cases considered. The lumped parameter simplification, by neglecting 
capacitance in the feed line, might cause the early time voltage cof Figs. 
3-4 and 3.5 to be slightly too large. And rinite conductivity would 
Slightly reduce the computed currents and voltages. 
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4.2 Discussion of Results 


LF vertical radiators respond quite dramatically to EMP from high 
altitude detonations. Peak currents calculated exceed 20 kA (Fig. 3.2) 
and peak load voltage is almost 5 MV (Fig- 3-5). Peak rate of rise is 
500 kA/usec (Fig. 3.1) and 30 MV/usec (Fig- 3-4). These values can be 
compared with average peak lightning currents of 10-20 kA and rate of 
rise of 5-20 kA/ysec. The duration of current surges from lighcning 
strokes would generally be longer than those induced by EMP. For these 
calculations we have not included any lightning arresters or other pro- 
tective devices. Such devices would shorten the duration of EMP induced 
surges. 

The effect of the angle of incidence 9 is quite pronounced. The 
early time current and voltage at 30° is almost twice what it is at 
90°, confirming the sin=! 8 behavior. On the other hand, end effects 
occur sooner at 30°, reducing the low frequency content. Thus the late 
time behavior is greater for 90°. In general it will be true for long 
antennas that effects cepending on initial rate of rise will be maxi- 
mized for small 6, whereas effects that depend on total surge energy 
are maximized for broadside incidence. 
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